function output = gammaup(gamma)

p0 = -2;

output = 2*abs(p0)*gamma - 9.8 - gamma*sqrt(2*abs(p0)*gamma);

end
